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The low energy excitations in Na2lr03 have been investigated using resonant inelastic x-ray scat- 
tering (RIXS). A magnetic excitation branch can be resolved, whose dispersion reaches a maximum 
energy of about 35 meV at the F-point. The momentum dependence of the excitation energy is 
much larger along the F — X direction compared to that along the V — Y direction. The observed 
dispersion relation is consistent with a recent theoretical prediction based on Heisenberg-Kitaev 
model. At high temperatures, we find large contributions from lattice vibrational modes to our 
RIXS spectra, suggesting that a strong electron-lattice coupling is present in Na2lr03. 

PACS numbers: 75.10.Jm, 75.30.Ds, 78.70.Ck 



The physics of iridates has drawn considerable atten- 
tion recently ii"— One of the reasons for this surge of in- 
terest is the fact that the local magnetic moment arises 
from a spin-orbit coupled jcff=l/2 state rather than a 
spin only state with quenched orbital moment as usually 
found in 3d transition metal compounds. One of the con- 
sequences of the Jeff =1/2 ground state is that the mag- 
netic interactions between such jcS moments can take on 
a form that is different from usual Heisenberg superex- 
change interactions. Specifically, bond-dependent Kitaev 
interactions are believed to arise when these jeff nio- 
ments reside on a honeycomb lattice. Due to the bond- 
dependence of the Kitaev interaction, strong frustration 
exists within this model which can induce a spin-liquid 
ground statei^i^ 

Na2lr03 is a promising candidate in which Kitaev in- 
teractions might be realized. In Na2lr03, edge-sharing 
IrOg octahedra form a honeycomb net which is deco- 
rated by Ir''+ ions with a 5d^ electronic configuration.*^ 
Experiments have indicated that Na2lr03 is a Mott in- 
sulator with an optical gap of '--^350 meVj^ Below T — 
15 K it orders antiferromagnetically in a so-called zig- 
zag structure!^ Recent resonant inelastic x-ray scatter- 
ing (RIXS) results have found evidence that the large 
spin-orbit coupling (SOC) in Na2lr03 is a dominant en- 
ergy scale.— This would cause the Ir moments to ac- 
quire a significant orbital component, giving it a jeff =1/2 
ground state,— and rendering the Kitaev interaction^'-*'' 
relevant. 

In particular, a combination of magnetic susceptibil- 
ity measurements and theoretical calculations on both 
Na2lr03 and Li2lr03 was used to claim that the spin 
model of this system is described as a Heisenberg-Kitaev 
model.— However, strong second and third neighbour 
exchange interactions are required to account for the 
ordered state, the large frustration parameter, and the 



magnon excitations at low temperatures.*^ Indeed, recent 
first-principle calculations suggest a considerable elec- 
tron delocalization in quasi-molecular orbits (QM0) ,*^i^° 
which is compatible with the experimental results. At 
present it is thus unclear, whether a localized scenario 
with contribution of Kitaev interactions or alternatively 
a more itinerant QMO picture is more appropriate for 
Na2lr03. One of the methods to test the existence of Ki- 
taev interactions in Na2lr03 is to study its magnetic ex- 
citation spectrum. RIXS at the Ir L3 edge has been suc- 
cessfully used to map "magnon" dispersions in iridates. 
For example, a sizable departure from a Heisenberg-like 
model was observed going from ST2lrO^ to Sr3lr207~ 
and was attributed to new bond-dependent magnetic in- 
teractions. 

Here we present Ir L3 edge RIXS results on a single 
crystal sample of Na2lr03. Our high- resolution RIXS 
measurements allow us to resolve new low energy excita- 
tions. In particular, at the F-point we observe a magnetic 
excitation which is centered around 35 meV. This excita- 
tion is almost dispersionless along the Q = [0 1 0] direc- 
tion but softens along Q = [1 0]. The large magnetic 
intensity observed at the F-point is consistent with the 
theoretical prediction including a bond-dependent Kitaev 
interaction in Na2lr03 4 An unusually strong electron- 
lattice coupling is also evident in Na2lr03 manifesting 
itself in the appearance of resonantly enhanced vibra- 
tional modes at high temperatures. 

The RIXS experiment was carried out at the Advanced 
Photon Source using the 30ID MERIX spectrometer. A 
spherical (2 m radius) diced Si(844) analyzer was used. 
The overall energy resolution (FWHM) in this configu- 
ration was ^ 30 meV. In order to minimize the elastic 
background intensity, most measurements were carried 
out in a horizontal scattering geometry near Q = (0 
6.75), for which the scattering angle 20 was close to 90°. 
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FIG. 1: (Color online) a) Wide range RIXS spectrum for a single- crystal sample of Na2lr03 at Q = (0 6.7) obtained with Ei 
= 11.217 keV. b) Detailed view of the low-lying RIXS excitations in a) taken at two different momentum transfers along the 
L-direction. Dashed lines are a fit to the elastic line (details provided in main text), c) and d) Incident energy dependence of 
the RIXS spectum. The intensity scale for (d) is shown. 



A single crystal of Na2lr03 was grown by the solid state 
synthesis method, as described in Ref. [Tj. The Na2lr03 
crystal was plate-like with a flat shiny surface, and a 
surface normal in the (001) direction. Throughout this 
paper we will use the C2/m notationi^i^i to describe the 
crystal structure. 

The RIXS process at the La-edge of Ir (or any other d 
electron system) is a second order process consisting of 
two dipole transitions {2p — >■ 5d followed by 5d -^ 2p) . As 
such, it can probe excitations between the d- levels ,-i^^"— 
collective magnetic excitations^iiSi2^i2S and even lattice 
vibrational modesi^"— 

In Fig. [H (a), a representative high-resolution RIXS 
spectrum of Na2lr03 is plotted on a wide energy scale. 
This scan was obtained at T = 10 K with Q = (0 6.7) 
and plotted as a function of energy loss {uj = Ei — Ef). 
The incident energy, Ei = 11.217 keV, was chosen to 
maximize the resonant enhancement of the spectral fea- 
tures of interest below 1 eV. Multiple peaks are observed 
between 0.4-1 eV (labeled A, B, and C). According to 
Ref. [la, features B and C correspond to excitations be- 
tween the jeff=3/2 and 1/2 states, and feature A corre- 
sponds to an exciton formed by a particle-hole pair across 
the gapji^ These features will not be discussed further in 
this paper. Below these peaks (w < 400 meV) we can 
observe the onset of the gap followed by what appears to 
be a large elastic line (w = 0). The elastic line is plotted 
on an expanded energy scale in Fig. [T] (b); an asymmet- 
ric line shape is evident with a maximum intensity at an 
energy loss of about 35 meV. Moving to a different mo- 
mentum transfer, Q=(0 6.5), the asymmetry remains, 
but since 29 is now further away from 90° the intensity of 
the elastic line has increased. By comparing the spectra 
at these two Q-points it is clear that the asymmetry is 
caused by new low- lying excitations, labelled Ag. To em- 
phasize this, the elastic contributions to the spectra are 
shown as dashed lines using the instrumental resolution 
function (described later in the text). 

Insights on the origin of Aq can be obtained through 



the incident energy dependence, which reveals which un- 
occupied 5d states comprise the intermediate state of the 
RIXS process. In Fig. [l](c) we compare spectra taken at 
Ei = 11.217 keV and E, = 11.220 keV. At Ei = 11.220 
keV, corresponding to exciting an Ir 2p3 /2 core electron 
into the unoccupied Ir 5d Cg level;^ we notice a drastic 
drop in Aq intensity. In Fig. [T] (d) we plot the inci- 
dent energy dependence of the RIXS spectrum obtained 
with a low resolution setup (FWHM ^ 150 meV). From 
the intensity plot we observe that Aq resonates around 
Ei = 11.217 keV, corresponding to the t2g intermedi- 
ate state, just like features at higher energiesJ^ In other 
words, the RIXS process for Aq, as well as those for A-C, 
involve an intermediate state which excites a 2p3/2 elec- 
tron into the 5d t2g states. The increased intensity of 
the elastic line below Ei — 11.211 keV comes from the 
decreased absorption as the incident energy falls below 
the Ir Ls edge. 

Let us consider possible explanations for the origin of 
Aq. We can quickly discard excitations between j^g lev- 
els since peak B represents the lowest energy transition 
possibleJ^ In addition, the peak position of Aq is an order 
of magnitude smaller than the optical gap^^ which ex- 
cludes any charge related excitations. This leaves us with 
either magnetic or lattice excitations,"^^ both of which are 
expected to appear in this energy range. We will model 
the total intensity of our RIXS spectrum with 3 compo- 
nents: 



/ = /''s + /' + /" 



(1) 



where I^^, P and 7™ represents the elastic background, 
lattice vibrational modes and magnetic excitations, re- 
spectively. We use a pseudo-Voigt lineshape as the in- 
strumental resolution function, which is a mixture of 
Lorcntzian and Gaussian functions with equal widths and 
amplitudes: 
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Here the FWHM was kept at the resolution-hmited value 
of 2r — 33 meV. For the rest of this paper we will refer 
to R{luo,uj) as the resolution function which is centered 
at ujQ. An example of this function is shown in Fig. [T] (b) 
and (c), where the elastic line has been fit using I''^{lu) — 
A^'' ■ R{0,uj), with A^' the amplitude of the elastic line. 
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FIG. 2: (Color online) a) Our fit function for the T = 400 
K RIXS spectrum, showing the number of resolution limited 
functions which represent the elastic line and the lattice vi- 
brational modes, b) Calculated total intensity of the Stokes 
(u > meV) and anti-Stokes {uj < meV) contribution to 
the lattice vibrations as a function of temperature, c) Tem- 
perature dependence of the RIXS spectrum taken at 10 K 
< T < 400 K for Q = (0 6.7). Spectra have been shifted 
vertically for clarity. The thick black solid lines show the elas- 
tic component, while the thin black lines include the lattice 
vibrational contribution as well. 

Since it is difficult to distinguish between /' aird I™, es- 
pecially when lattice and magnetic excitations appear on 
the same energy scaled we rely on temperature depen- 
dence to analyze the observed spectra. At temperatures 
well above the characteristic magnetic energy scale, lat- 
tice excitations will dominate the inelastic signal, while 
at lower temperatures the vibrational modes will coexist 
with spin excitations. In Fig. [5] (b) we show the RIXS 
spectrum taken at T=400 K, well above the Curie- Weiss 
temperature, 0cw — -116 K, of Na2lr03J^ At this el- 
evated temperature only lattice vibrations are expected 
to contribute to the low energy RIXS signal. Note that 
this high temperature spectrum is much more symmetric 
than the T = 10 K data. The maximum intensity posi- 
tion seems to be closer to the elastic line due to the large 
spectral weight seen on the energy gain side. 

In a recent RIXS study of the edge sharing cuprate 
compound Ca2+5a;Y2-5xCu50io,— it was found that the 
resonant enhancement of vibrational modes does not oc- 
cur uniformly. That is, certain modes and their higher 
harmonics are selectively enhanced in a RIXS experi- 



ment. Motivated by this work we model the I contri- 
bution with a series of vibrational modes: 



l'{u) = ^ 4 MuJ,) + 1)R{lu„lu) + n{u,)R{-u,, u)] , 

(3) 
where the i-th harmonic of the vibrational mode is de- 
scribed by a resolution limited peak centered at Wi. 
n{u)i) = \j[e^il^BT _ j^-j jg |-j^g Bose population factor. 

Note that the principle of detailed balance constraints 
the ratio of Stokes and anti-Stokes peaks, leaving the 
amplitudes (A') as only adjustable parameters. Based 
on recent optical conductivity data'^A we fixed our lowest 
mode at wi — 18 meV— with higher levels placed in a 
harmonic orderi^i^^i^l of the first peak (uj^ — 36 meV, 
W3 == 54 meV and uj/^ = 72 meV) . In Fig. [2] (a) we show 
our fit result as a solid black line using / = I^^ +P, where 
only A'^'' and A\ were fit. 

Having estimated / = I^^+P at T = 400 K we can pro- 
ceed to calculate its expected temperature dependence 
for T<400 K by assuming that I^^ is independent of tem- 
perature and that the intensity of the vibrational modes 
follows the thermal population factors: (^(0;^) -f 1) and 
n{uJi) for Stokes and anti-Stokes, respectively. One can 
see the expected temperature dependence of the vibra- 
tional modes intensity in Fig. [11(b). At T = 10 K the 
anti-Stokes vibrational mode contribution vanishes. In 
Fig. [5] (c) the calculated intensity of / = I^s + p is plot- 
ted as thin solid lines. The RIXS spectra shows extra 
intensity for T< 200 K. This extra intensity grows as 
temperature decreases and leads to the appearance of Aq 
(peaked around 35 meV). 

This study of temperature dependence clearly shows 
that on top of the vibrational modes, there exists addi- 
tional intensity that grows with decreasing temperature. 
The most likely origin of this intensity is magnetic for 
the following two reasons: 1) the temperature scale for 
the onset of this intensity is the same order of magni- 
tude as 6 
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and 2) the ground state of this sample 
is magnetically orderedJ^ We also note that these lat- 
tice vibrational modes are unusually strong, with higher- 
harmonics carrying appreciable spectral weight which is 
even comparable to the charge excitations seen in Fig. 
[1] (a). In Ir L3 edge RIXS experiments electron-lattice 
coupling can cause lattice vibrational excitations to ac- 
quire significant spectral weight;— This comes about due 
to the sudden change in the charge density around the Ir 
atoin during the 2p — J> bd absorption process. The large 
inteirsity of the vibrational modes and significant spectral 
weight of higher-harmonics in Na2lr03 is thus an indica- 
tor of strong electron-lattice couplingj^i2? which might 
help explain the unusual spectral broadening observed by 
angle resolved photoemission spectroscopy.— 1^ 

We now move on to discuss the dispersion of the mag- 
netic excitation. In Fig. [3] (a), the momentum depen- 
dence of Aq along both the T — X and T — Y directions 
(from (0 0) to (0.5 0)/(0 0.5)) at T = 10 K is shown, with 
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FIG. 3: (Color online) (a) Momentum dependence of the 
RIXS spectra taken along the T — X and T — Y directions at T 
= 10 K. Superimposed on each spectra is the background con- 
tribution from the elastic line and the vibrational modes, (b) 
The magnetic signal after subtracting the background, and 

(c) a direct comparison between Q = (0 0) and Q = (0.5 0). 
In (a) and (b) spectra have been shifted vertically for clarity. 

(d) Fit results for the peak position of the magnetic signal in 
(b). 



spectra shifted vertically for clarity. Also shown in this 
figure as dashed lines are the contributions from (I^^+P) 
in order to show the magnetic contribution clearly. We 
note that the intensity of the elastic line was allowed to 
vary between different values of Q, which is expected in 
general (for a diffuse scattering intensity). The tail on 
the energy gain allows us to estimate the elastic inten- 
sity. This is possible because the inelastic (lattice vibra- 
tions) contribution to the energy gain side at this tem- 
perature is negligible. However, the contributions from 
the vibrational modes were fixed. At both Q = (0.38 
0) and Q = (0.5 0) the spectral weight of Aq shifts to 
significantly lower energies. In Fig. [3Kb), the magnetic 
intensity has been plotted by subtracting I^^ + P contri- 
butions from the raw spectra. The large error bars on the 
energy gain side reflect the uncertainty arising from the 
elastic background fitting. Despite this, it is clear that 
at the F-point Aq forms a well defined feature which is 
centered around 35 meV, and that Aq disperses towards 



lower energy along the Q = [1 0] direction. At Q = (0.5 
0) most of the spectral weight from the magnetic mode 
is only visible below 30 meV. Fig. [3] (c) shows the differ- 
ence spectrum obtained at Q = (0 0) and Q = (0.5 0) 
without an offset. Although the large error bars make 
it difficult to extract the shift of the peak, a decrease of 
roughly 25 meV in energy is observed. In an attempt to 
capture this dispersion we fit the magnetic signal using 
the resolution function provided in Eq. (2). The fitted 
peak positions are plotted in Fig. EJd), showing signif- 
icant dispersion from T to X . On the other hand, the 
momentum dependence along the T — Y direction (from 
Q = (0 0) to (0 0.5)) is much weaker (see Fig. [3](b) bot- 
tom spectra); no significant shift is observed. This is in 
stark contrast to the 25 meV dispersion observed along 
the T — X direction. 

Our observation of a high energy ('^35 meV) mag- 
netic excitation is rather surprising. In a recent inelas- 
tic neutron scattering experiment on a powder sample of 
Na2lr03 a magnon mode below 6 meV was identifiedj^^ 
A pure Heisenberg model with antiferromagnetic inter- 
actions and additional long range exchanges was found 
to be adequate in describing this result J^ The calculated 
dispersion of the high energy branch at the F-point, how- 
ever, only reaches about 5 meV (see supplemental ma- 
terial in Ref. llSh . which is significantly lower than the 
position of Aq. Given the large energy separation, it is 
difficult to explain both sets of data with a purely Heisen- 
berg Hamiltonian. Recently, Chaloupka et al.*^ were able 
to explain the observed neutron data^^ by adding a Ki- 
taev term to the Heisenberg Hamiltonian. Moreover, 
this additional Kitaev term would generate a high energy 
magnon branch. According to Ref. |8), this branch would 
reach 20 meV at the F-point with anisotropic dispersion 
along the T — X and T — Y directions. These predictions 
therefore seem quite consistent with our experimental ob- 
servations in Na2lr03. We acknowledge that in order to 
extract the size of the Kitaev term, and more impor- 
tantly to determine its sign,^'^® higher resolution RIXS 
data will be required. 

In conclusion, we have identified new low-lying excita- 
tions in Na2lr03 using high-resolution Ir L3 edge reso- 
nant inelastic x-ray scattering. Temperature dependence 
reveals two distinct modes: a dominant lattice vibra- 
tional mode at high-temperature and a magnetic exci- 
tation which appears below T = 200 K and reaches max- 
imum intensity at T = 10 K. The vibrational excitations 
were fit using an optical phonon mode at wi — 18 meV, 
with noticable spectral weight on the next three higher- 
harmonics. This suggests that the electron-lattice cou- 
pling is very strong in Na2lr03. The magnetic mode 
shows peculiar momentum dependence at T = 10 K, 
reaching a maximum energy of ~ 35 meV at the F-point 
and dispersing to lower energies along the T—X direction. 
The observed dispersion of magnetic excitation is consis- 
tent with theoretical calculations based on a local spin 
model with both Heisenberg and Kitaev interactions. 
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